Entanglement by a beam splitter: nonclassicality as a prerequisite for entanglement 
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A beam splitter is a simple, readily available device which can act to entangle the output optical 
fields. We show that a necessary condition for the fields at the output of the beam splitter to be 
entangled is that the pure input states exhibit nonclassical behavior. We generalize this proof for 
arbitrary (pure or impure) Gaussian input states. Specifically, nonclassicality of the input Gaussian 
fields is a necessary condition for entanglement of the field modes with the help of the beam splitter. 
We conjecture that this is a general property of the beam splitter: Nonclassicality of the inputs is 
a necessary condition for entangling fields in the beam splitter. 

PACS number(s); 03.67.-a, 89.70.+C 



I. INTRODUCTION 

Entanglement is at the heart of current development 
of quantum information processing Q. Entanglement- 
assisted communication can enlarge the channel capacity 

and enhance channel efficiency Entanglement may 
play a key role in secure communication Q . In quantum 
computation, of course, qubits are massively entangled. 

The generation and characterization of entanglement 
has been studied extensively. In particular, a recent ex- 
perimental advance realized the generation and distilla- 
tion of polarization-entangled photons toward optimal 
entanglement in a 2 x 2-dimensional Hilbert space [^. 
The polarization-entangled photons are generated using 
type I or type II parametric down conversion. The para- 
metric down conversion is also a standard technique to 
produce a two-mode squeezed state, which is an entan- 
gled state in an infinite dimensional Hilbert space ||^. 

The beam splitter is also one of only a few 
experimentally-accessible devices which may act as an 
entangler. There have been some previous studies of a 
beam splitter as an entangler 1^-^. In particular, Paris 
studied entanglement properties of the output state 
from a Mach-Zehnder interferometer for squeezed input 
states. The action of a linear directional coupler can 
also be described by the beam splitter operator. Photon 
statistics and nonclassical properties of the output fields 
from a linear directional coupler were studied for Fock 
and squeezed inputs pOt . 

In this paper we investigate the entangling properties 
of a beam splitter for various pure input states including 
Fock states and squeezed states. We find a simple for- 
mula to determine the entanglement of output fields for 
squeezed input fields. We also study the entanglement 
of output fields when the input fields are in a Gaussian 
mixed state and provide a sufficient condition for input 
fields to have no entanglement in the output state: when 
two Gaussian "classical" fields are input to the beam 
splitter, the output state is never entangled. We find 



that entanglement of the output state is strongly related 
to the nonclassicality of the input fields. 



II. FOCK-STATE INPUT 

Fig. 1 shows the schematic arrangement of a beam 
splitter. The input field described by the operator a is 
superposed on the other input field with operator & by a 
lossless symmetric beam splitter, with amplitude reflec- 
tion and transmission coefficients r and t. The output- 
field annihilation operators are given by 



c = BaB^ , d = 
where the beam splitter operator is lO 



B = exp 



-(d^Se'*-a&^e- 
2^ 



(2.1) 



(2.2) 



with the amplitude reflection and transmission coeffi- 
cients 



t — cos - 
2 



r = sm — . 
2 



(2.3) 



The beam splitter gives the phase difference (f) between 
the reflected and transmitted fields. 

In this paper we are interested in entanglement prop- 
erties of the output state. Suppose that input states are 
two independent Fock states, |ni,n2) = \n\)a\'n2)h- The 
output fields are then a superposition of two-mode Fock 
states: 
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with 6 a Kronecker delta function. When the total num- 
ber of input photons is = ni + ^2, the output state 
becomes an [N + l)-dimensional entangled state. 

The von Neumann entropy is a measure of entan- 
glement for pure bipartite states (See e.g. [Q), which 
becomes ln(A^ -I- 1) when an {N + l)-dimensional bi- 
partite system is maximally entangled. The von Neu- 
mann entropy E{pa) for the reduced density operator 
Pa = TibB\ni,n2){ni,n2\B^ is 



E{Pa) 
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Fig. 2 shows the von Neumann entropy E{pa) is a 
function of the reflection coefficient r and configuration 
of input photon numbers. It is interesting to note that 
the entropy does not necessarily maximize for a 50:50 
beam splitter. This is discussed further in the following 
subsections. 



A. SU(2) coherent state 

When N number of photons are injected into one input 
port while no photon is injected into the other input port, 
the output state turns into a state generally known as an 
SU(2) coherent state Substituting m = and 



n2 ^ N into Eq.(2.4), we find the SU(2) coherent state 



N 



B\0,N)=J2c^!\k,N-k) 



where 



(2.7) 
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The von Neumann entropy for the reduced density op- 



erator p is E{pa) = Y.k=' 



k=0 l^fe 



In 



In Fig. 2, 



the von Neumann entropy for = 10 is plotted, which 
shows that the measure of entanglement is a convex func- 
tion with its maximum for a 50:50 beam splitter, i.e. 
r = t = 1/2. In particular, when = 1 the output state 
is l/\/2(|0, 1) + |1,0)) for a 50:50 beam splitter @,||. 



B |?i,7i) = ^ e 
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2m — y nl I \ / 

for a 50:50 beam splitter. This shows that the possibil- 
ity of having odd numbers of photons is zero |0). This 
is an extension of the well-known result of Sjl,!) = 
^(|0,2) -I- e*"^|2,0)) Output state |1,1) may re- 

sult from transmission or reflection of both the photons. 
The two cases destructively interfere to remove the |1, 1) 
state in the output state. In fact, the output state is the 
maximally entangled state in the Hilbert space composed 
of |0) and |2). We now see why entanglement does not 
maximize when the same number of photons are injected 
to a 50:50 beam splitter. This is due to the fact that 
odd-number states destructively interfere and do not ap- 
pear in the output state. It is also true that the output 
state can be considered in {n -I- l)-dimensional Hilbert 
space composed of |0), |2), • • •, \ 2n) instead of {2n -1-1)- 
dimensional space. 

With the use of a beam splitter, there are two ways to 
generate entangled states in (n -I- l)-dimensional Hilbert 
space. One way is to put total of n photons into a beam 
splitter and the other way is to put n photons into each 
input port of a 50:50 beam splitter. By comparing the 
von Neumann entropies for the both cases, we find that 
the latter case of using a 50:50 beam splitter does not 
bring about the best entanglement. 



III. SQUEEZED STATE INPUTS 

Generating Gaussian states, in particular, coherent 
states and squeezed states has become a standard experi- 
mental technique. When two coherent states are incident 
on a beam splitter, the output is given by 

Bba{a)Db{P)% 0) = Da{ta + re'^ l3)Db{tl3 - re-''^«)|0, 0) 

(3.1) 

where D{a) — exp(aat — a* a) i s the displacement opera- 
tor ||l^. The output state (B.l) is clearly not entangled. 
It is further found that displacing the input fields does 
not increase entanglement of the output fields because 
the impact of displacing the input fields can always be 
canceled by local unitary operations on the output fields. 

When the two input fields are squeezed, the output 
state from a beam splitter is 



B. Input fields of same number of photons 

In Fig. 2, it is interesting to note that, for a 50:50 beam 
splitter, the entanglement shows a dip when ni — n2. 
When the two input Fock states have the same number 
of photons, i.e., ni = n2 = n, the output state is 



i?^a(Cl)^6(C2)|0,0) 



where the squeezing operator |18 



^(C) = exp ( ic*«' - \Ca^^ 



(3.2) 
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with the complex squeezing parameter C = sexp(?(^). 
The phase of the squeezing parameter determines 
the direction of squeezing. Using the rotation operator 
— exp{i'da'' a) the following can be written 



5(0,0)5(0 = B{0,<p)R{ip/2)S{s)R\ip/2) 

= R{ip/2)B{e, 4> - ip/2)S{.s)R\ip/2), 



(3.4) 



where, in order to specify the parameters 9, (p of the 
beam splitter operator, the beam splitter operator has 
been denoted by B(9, cj)). The first rotation operator in 
the last line of Eq. ( |3.4| ) is canceled by local operation 
and the last rotation operator does not change the state 
when it applies to the vacuum. Now, we have found that 
the relative phase between the amplitude reflection and 
transmission coefficients gives the effect of the rotation of 
the squeezing angle for the input fields. Without losing 
generality, we take the input squeezing parameter to be 
real while keeping variable. 

Th e von Neumann entropy E{pa) of the output state 
(3.2) is plotted in Fig. 3 against squeezing parameter S2 
and reflection coefficient for si — 0.5. The relative phase 
(/) = in Fig. 3(a) and 7r/2 in Fig. 3(b). We find that the 
entanglement of the output state depends on the degrees 
of squeezing for input fields and the reflection coefficient. 
We also note that the relative phase <j) hence relative an- 
gle of squeezing for input fields plays an important role. 
For a 50:50 beam splitter, the entanglement of the out- 
put state is minimized when = 0, while it is maximized 
when (j) — n /2. In other words, for = 0, the entangle- 
ment of the output state is maximized if the two input 
fields are squeezed along the conjugate quadratures in 
phase space. To analyze the output state ( ^ ) further, 
consider the following relation for a 50:50 beam splitter 
of = £7r/2 (£ = 0,1,2,- • •). In this case, the output 



state (3^) can be written as 



B(V4,0)5a(si)5b(s2)|O,O) 

= 5aQ(si + S2e2^*))5fcQ(sie-2'^ 
xS„b(^i(sie'*-S2e-^*)j |0,0) 



+ S2) 



(3.5) 



where Sab[0 = exp(— ^afo -I- C,*a)w) is the two-mode 
squeezing operator. The single-mode sq ueez ing operators 
Sa and Sb in the right-hand side of Eq. ( |3.5| ) do not con- 
tribute toward entanglement of the output state because 
they can be canceled by local unitary operations. Thus 
only the two-mode squeezing operator Sab determines the 
entanglement of the output state as only it represents a 
joint action on both pairs of the bipartite system. For 
a given squeezing, si and S2, when = 7r/2, the output 
state is maximally entangled. When = 0, entanglement 
is minimized. In fact, if si = S2 we completely lose entan- 
glement for = 0. We notice that a two-mode squeezed 
state is produced from a single-mode squeezed state by an 
action of a beam splitter and local unitary operations. In 



contrast to the case of the Fock-state input, the relative 
phase between reflection and transmission plays an im- 
portant role for the case of squeezed input fields. 

So far, we have studied only pure input states. From 
what we have learned we can conclude that the nonclas- 
sical behavior of the input fields is a necessary condi- 
tion for the output fields to be entangled. Specifically, 
the only pure state which does not possess nonclassical 
properties is a coherent state (Its P-function is positive 
well-defined. See the discussion in the next section.). As 
it is well known coherent inputs never become entangled 
in the beam splitter, that is the output can always be 
written in the factorized form. On the other hand, as we 
have shown above, nonclassicality of the inputs is not a 
sufficient condition for the entanglement. 



IV. GAUSSIAN MIXED STATE INPUT 

When the input fields are mixed, the output fields 
from a beam splitter are also mixed. A general mixed 
continuous- variable state is not easy to deal with because 
of its complicated nature. However, for a Gaussian two- 
mode state, the separability condition has been studied 
extensively [|l^-|l)- 

The separability of a Gaussian state is discussed with 
quasi-probability functions and their characteristic func- 
tions in phase space. There are a group of quasi- 
probability functions including the Wigner function, the 
Husimi Q function and the P-function [p2| . In partic- 
ular, the P-function can be used as the measure of the 
nonclassicality of the given field. For example, if a single- 
mode state is nonclassical its density operator p'^' can be 
written as 



p"' = / P{a)\a){a\d'a 



(4.1) 



where the P-function P(a) is positive and well-behaved. 

It has been shown that if a two-mode Gaussian state 
is represented by a positive well-behaved P-function 
P{a,l3), the state is separable 21|. Suppose two clas- 
sical states of P-functions Pa{ a) a nd Pb{f3) are incident 
on a beam splitter. Using Eq.(4T), the density operator 
for the output state is written as 

B J Paia)Pbma)a{a\®\(3)a{f3\d^ad^pB^ 

= J Paia)Pbil3)\ta + re"^l3)a{ta + re''^l3\ 

(g) I - re-^^a + tf3)b{-re-">'a + tl3\d^ad^l3 
= I Pa{tl-re'H)Pb{re-'^^ + t6Ma{l\®\5)a{5\d^ld^5. (4.2) 



Here Pa(t7 - re''^5)P6(re"'"^7 + t5) is the two-mode P- 
function for the output state. Because Pa (a) and Pb{/3) 
are positive well-defined under the assumption of classical 
input fields, Pa(i7 — re^'^S)Pb{re'^'^j -\- tS) is also positive 
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well-defined. We have proved a sufficient condition for 
separability of the output state from a beam splitter: 
when two classical Gaussian input fields are incident on 
a beam splitter, the output state is always separable. It 
follows that for creating a Gaussian entangled state with 
a help of a beam splitter, it is necessary that the input 
exhibits nonclassical behavior. 

We have already seen that two nonclassical input fields 
do not necessarily bring about entanglement in the out 
state as two squeezed state inputs may not entangled in 
the beam splitter. We investigate the entanglement of 
the output state when two Gaussian mixed states are 
incident on a beam splitter. 

The necessary and sufficient criterion for the separa- 
bility of a Gaussian mixed state has been studied us- 
ing the Weyl characteristic function C^"''(C,?y) [^^-^. 
For a two-mode Gaussian state of density operator 
Pab, the Weyl characteristic function |^^, C'''^\(,r]) = 
TipabDaiODbiTj), can be written as 



C('^)(C,r;) =cxp 



(4.3) 



where M is a 4 x 4 matrix which completely determines 
the statistical properties of the Gaussian state. Duan et 
al. [|2^ found that after some local operations, it is possi- 
ble to transform the state into another that is represented 
by the matrix, 



M' 



(hi ci 

&2 C2 

ci di 

V C2 d2, 



(4.4) 



where parameters &i, di and Ci satisfy 

hi - 1 ^ &2- 1 
di - 1 ^ - 1 



|C1|-|C2| = ^{hi-\){di-\) 

-\/(&2 - l)(d2 - 1). 



(4.5) 



(4.6) 



Note that parameters ci^2 determine the correlation be- 
tween two modes. The necessary and sufficient criterion 
for separability reads then 



{{Auf) + {{Mf) > ql 



1 



(4.7) 



where — \J [di — l)/(fci — 1) and two operators u and 
V are defined as 



—;={a + a') - - — - ^ 
V2 \ci\V2qo 



iqo 
V2 



(at 



C2 



|C2| ^/2qo 



(6t + h) 



(4.8) 



When two mixed states of density operators pa and 
Pb are input to a beam splitter, the density operator for 



the two-mode output field is pout = BpW . The Weyl 
characteristic function for the output field is 

C'ZI (C, V) = (iC + reS)C\r^ i-re-^^C + tv) (4.9) 

which is obtained using the relation W Da{C)Db{ri)B = 
DaitC + re"t'7])Db{-re-''t'C + tr]). 



A. Squeezed thermal state inputs 

Consider two thermal states of the same average pho- 
ton number n. The density operator for the thermal field 



Pth 



^ (»)" , w I 
> -, zv-. — \n){n\. 



(4.10) 



Suppose the thermal fields are respectively squeezed be- 
fore they are mixed at a beam splitter. From the earlier 
section, we know that two squeezed vacua result in max- 
imum entanglement for the output field when </> = 7r/2. 
We thus restrict our discussion to the case cj) = it/2 for 
the study of two squeezed thermal state inputs. We also 
assume that the incident fields are equally squeezed. 

The squeezed thermal field S{s)pthS^ (s) is represented 
by the following characteristic function: 



exp 



-(2n- 
2^ 



l)e-^Xf 



(4.11) 



The squeezed thermal state is said to be nonclassical 
when one of the quadrature variables has its variance 
smaller than the vacuum limit; the squeezed thermal 
state of (4.11) is nonclassical when [p5[ 



{2n 



Dc-^" - 1< 0. 



(4.12) 



Throughout the paper s > is assumed without loss of 
generality. 

For the maximum entanglement of the squeezed input, 
let us consider a 50:50 beam splitter. Substituting cj^J 
of (4.11) into Eq.(|4.9|), the matrix elements in Eq.(| 



are found: 



bi 



Cl 



C2 



di 

1 

2 
1 

2 



1 



{2n- 
{2n- 



2 

l)(e-2^ 



(2n+l)(e^^ + e-^^) 



l)(e^ 



2.S 



e^^) 



(4.13) 



The separability condition (4.7) in this case reads that 
the output state is se parab le when &i — 1 > |ci|. Sub- 
stituting bi and ci in ( 4.13| ), it is found that the output 
state i s sep arable when {2n + l)e^^'' — 1 > 0. With help 
of Eq.(4.12), we write that the output state is entangled 
when the squeezed thermal input fields becomes nonclas- 
sical. 
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B. Squeezed thermal and vacuum input states 

Suppose a squeezed thermal state is incident on one 
input port and vacuum is incident on the other input 
port. As was done earher, we assume (/) — 7t/2 for the 
beam sphtter. In this subsection we release the condition 
of the 50:50 beam splitter hence the output state depends 
on the reflection coefficient of the beam splitter. The 
output state is then represented by the matrix M with 
its elements: 



bi = {2n + l)r^ + e-^H^ 



-2s 



bi = r^(2n+ f)e 
di = f{2fi + l)c-'^' - 
ci = tr[{2n + f )e"2" 



62 = r^2n 
C2 — tr[{2n 



f )e2^ + 



- fl 



(4.14) 



The separability criterion (4.7) takes different forms de- 
pending on the positivity of 61 — 1 and di — 1 due to the 
definition of qo- When 5i — f > and c?i — f > 0, the 
separability criterion becomes 



V/(61 - l)(dl - f) + ^{bl - f)(di - f) > |ci| + |C2|. 

(4.15) 

Otherwise the separability criterion is 

- - - 1) + ~ - 1) > |C1| + |C2| 

(4.16) 



With the use of 61 and di in (4.14), we find that both con- 
ditions 61 - 1 > and di - 1 > i mply (2n-f l)e-2« - 1 > 
0. In this case, the inequality ( 4.15 ) is always satis- 
fied and the output state is separable. Howeve r, wh en 
{2fi + l)e^^'' — 1 < 0, the separability criterion ( 4.16| ) is 
never satisfied and the output state is entangled. Here, 
we confirm our earlier finding that the nonclassicality of 
the input state provides the entanglement criterion for 
the output state. When a squeezed thermal state and 
vacuum are incident on a beam splitter, the output state 
is entangled only if the squeezed thermal state is non- 
classical. 



C. Squeezed vacuum and thermal input states 

So far, we found that nonclassicality of the incident 
field plays an important role in the entanglement of the 
output field. Let us suppose that an input field is a 
squeezed vacuum and the other input field is a thermal 
state. Differently from the earlier cases in this section, 
one of the input states is always nonclassical while the 
other is always classical. Substituting the characteristic 
func tion s for the thermal state and squeezed state into 
Eq.(4.9), the char acte ristic function for the output field 
is represented by ( |4.3| ) with matrix M in the form (4.4) 
and the matrix elements are 



di = {2n 
Ci = tr{2n 



I)t2 + e-2v2 



62 = (2n -f l)r' + e'^H'^ 
d2 = (2n-f I)t2+e2v2 

C2 =tr(2n-Kl-e^^). (4.17) 



These elements do not satisfy conditions (4.5) and ([l.q). 
In order to use the separability criterion (4.7), the output 
state has to be locally transformed. 

Suppose the output fields are squeezed locally. Assum- 
ing equal degree of squeezing, s, for each mode, the trans- 
formed state is represented by ptrans = SaSbPoutSlSl- 
We use the identity S"^ {s)D{a)S{s) = D{are^ + ia^e""*), 
where sub-indeces r and i respe ctive ly denote real and 
imaginary parts, and definition (4.3), to find the Weyl 
characteristic function for the transformed state: 



where C^^'j is the characteristic function for the out- 
put state. After a little algebra, we find that the ma- 
trix elements representing C'(J^''„g(C, ry) is the same as 
those in (4.14) for the output state from a beam splitter 



Ciui iCre' + ^Ge-^ r?,e^ + zry^e"^) (4.18) 



when the squeezed thermal and vacuum are input fields 
but with squeezing factor —s. The separability criterion 
{2h + 1)6^^^* — 1 < 0, thus, applies for the output state 
when the two input fields are the squeezed vacuum and 
thermal field. The separability criterion coincides with 
the nonclassicality condition for the output field of mode 
c in Fig. 1. 



V. REMARKS 

We have considered the nature of the entanglement of 
output fields from a beam splitter for pure state inputs 
and for mixed Gaussian state inputs. In the case of pure 
states we have found that for Fock state inputs, the beam 
splitter is a tool to produce a (A^ -I- l)-dimensional entan- 
gled state, where N is the total excitation of the input 
fields. For squeezed vacuum inputs, the entanglement of 
output fields depends on many factors including the rel- 
ative angle of squeezing between two input fields. When 
the relative angle is appropriately chosen, the entangle- 
ment of the output state is maximized for a 50:50 beam 
splitter. From these results it directly follows that non- 
classicality of input pure states is a necessary condition 
for having entangled states at the output of the beam 
splitter. 

In the case of mixed states the analysis is more com- 
plicated since there does not exist a necessary and suf- 
ficient condition for inseparability of arbitrary infinite- 
dimensional bi-partite systems. Since the condition ex- 
ists for Gaussian states, we have concentrated our atten- 
tion on these states. We have proved a sufficient condi- 
tion for the output state of a beam splitter to be separa- 
ble (that is they are not entangled): if both the Gaussian 
input fields are classical, it is not possible to create entan- 
glement in the output of the beam splitter. From here it 
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automatically follows that nonclassicality is a necessary 
condition for the entanglement. 

These observations make us conjecture that nonclas- 
sicality of at least one of the input fields is a necessary 
condition for the output to be entangled. That is the non- 
classicality of individual inputs can be traded for quan- 
tum entanglement of the output of the beam splitter. 
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FIG. 1. Configuration of beam splitter operation. 
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(a) 




FIG. 2. The measure of entanglement E{p) is plotted 

using the von Neumann entropy for the reduced density oper- 
ator of the beam-splitter output field. The Fock-state input 
fields |fc, N — k) have total photon number N = 10. R = r^. 
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FIG. 3. The measure of entanglement E{p) for the 
beam-sphtter output field is plotted using the von Neumann 

entropy for the reduced density operator of the output field. 
The squeezing parameter for one squeezed input is fixed to 
si = 0.5 while the squeezing parameter for the other squeezed 
state is varied from S2 = to 1. The transmittivity R. The 
beam splitter gives phase difference = (a) and (j) = ■k/2 
(b) between the reflected and transmitted fields. R = r^. 
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